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Fluctuations of the chemical composition of the hadronic system produced in nuclear collisions are 
discussed using the $— measure which has been earlier applied to study the transverse momentum 
fluctuations. The measure is expressed through the moments of the multiplicity distribution and 
then the properties of are discussed within a few models of multiparticle production. A special 
On : attention is paid to the fluctuations in the equilibrium ideal quantum gas. The system of kaons and 

pions, which is particularly interesting from the experimental point of view, is discussed in detail. 
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' There are many sorts of hadrons produced in high energy collisions. The ratios of multiplicities of particles of given 
species to the total particle number characterize the chemical composition of the collision final state. The composition 
is expected to reflect the collisions dynamics. Generation of the quark-gluon plasma in heavy-ion collisions was argued 
long ago jlj to enhance the strange particle production. While the significant strangeness yield enhancement has been 
experimentally observed in the central nucleus-nucleus collisions at CERN SPS (see the data compilation and the 
recent review Q), it is a matter of hot debate whether the observation can be indeed treated as a plasma signal. 
The strange hadron abundance is naturally described within the models assuming the plasma occurrence at the early 
collision stage (see e.g. but the models, which neglect such a possibility (see e.g. HQ or the review §), can be 

On ■ also tuned to agree with the experimental data. Thus, it would be desirable to go beyond the average particle numbers 
and see whether the strangeness enhancement in the central heavy-ion collisions is accompanied with the qualitative 
change of the strangeness yield fluctuations. The equilibrium quark-gluon plasma scenario is obviously expected to 
lead to the smaller fluctuations than the nonequilibrium cascade-like hadron models but the specific calculations 
are needed to quantify such a prediction. Anyhow, it seems to be really interesting to study the strangeness yield 
fluctuations on the event-by-event basis. However, we immediately face the difficulty how to quantitatively measure 
the fluctuations in the events of very different multiplicity. The problem appears to be of more general nature. 

There are several interesting proposals to use fluctuations as a potential source of valuable information on the 
collision dynamics. If the hadronic system produced in the collision is in the thermodynamical equilibrium, the 
temperature and multiplicity fluctuations have been argued to determine, respectively, the heat capacity p|,p!o[ and 
compressibility ]il[ | of the hadronic matter at freeze-out. An extensive discussion of the equilibrium fluctuations can be 
found in In the experimental realization of such ideas one has to disentangle however the 'dynamical' fluctuations 
of interest from the 'trivial' geometrical ones due to the impact parameter variation. The latter fluctuations are 
very sizable and dominate the fluctuations of all extensive event characteristics such as multiplicity or transverse 
energy. The variation of the impact parameter can also influence the fluctuations of the intensive quantities e.g. the 
temperature. 

A specific solution to the problem was given in our paper fDjfl , where we introduced the measure of fluctuations or 
correlations which has been later called <&. ft is constructed in such a way that 4> is exactly the same for nucleon- 
nucleon (N-N) and nucelus-nucelus (A-A) collisions if the A-A collision is a simple superposition of N-N interactions. 
On the other hand, $ equals zero when the correlations are absent in the collision final state. The method proposed 
in |l3f| has been recently applied to the NA49 experimental data. The fluctuations of transverse momentum found in 
the central Pb-Pb collisions at 158 GeV per nucleon have appeared to be surprisingly small |L4|,[L5). ft has been also 
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claimed 15 that the correlations, which are of short range in the momentum space, are responsible for the nonzero 
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positive value of <& PT being observed. Our calculations of $ PT in the equilibrium ideal gas show |16| that <J> PT is positive 
for bosons, negative for fermions and zero for classical particles. When the hadronic system at freeze-out is identified 
with the pion gas, the calculated $ PT slightly overestimates the experimental value [pl| but the inclusion of the pions 
which come from the resonance decays removes the discrepency. An interesting analysis of the py— fluctuations within 
the so-called non-extensive statistics is given in |l7j. 

The theoretical analysis of the result [fLtjjlq] has provided a new insight into the collision dynamics. It has been 



argued |18| within the UrQMD model that the secondary scatterings are responsible for the dramatic correlation 
loss in the central collisions of heavy-ions. This conclusion however seems to contradict the results of the analysis 
Jl9| where the LUCIAE event generator has been used and the rescatterings are shown to reduce insignificantly the 
Pt~ correlations measured by <&. While the effect of the secondary interactions needs to be clarified, the smallness of 
the fluctuations observed in the central heavy-ion collisions [ pTp^ ] is a very restrictive test of the collision models. 
The so-called random walk model is ruled out because it gives much stronger correlations in A-A than N-N case p2| . 
The same holds for the LUCIAE event generator when the jet production and/or the string clustering is taken 
into account even at a rather moderate rate. On the other hand, the quark-gluon string model seems to pass the test 
successfully 

As argued in p2|| , the measure $ can be also applied to study the fluctuations of chemical composition of the 
hadronic system produced in the nuclear collisions. The chemical fluctuations seem to be even more interesting than 
those of the kinematical variables such as px- The final state momentum distribution of hadrons characterizes the 
system at the moment of freeze-out, while the system chemical composition is fixed at the earlier evolution stage 
- the chemical freeze-out when the secondary inelastic interactions are no longer effective. The total strangeness 
yield presumably saturates even earlier and the subsequent interactions are mostly responsible for the strangeness 
redistribution among hadron species. 

The NA49 Collaboration plans to study the chemical fluctuations in heavy-ion collisions at CERN SPS (23). Since 
the <f>— measure will be used in these studies, it is desirable to better understand the properties of <& when applied to 
the chemical fluctuations. This is the aim of our note. At the beginning we express $ through the commonly used 
moments of the multiplicity distribution and analyze the result within several models of the distribution. Then, we 
compute the measure for case of the two-component ideal quantum gas in equilibrium. The system of kaons and 
pions is discussed in detail. In particular, the role of resonances is analysed. 

Let us first introduce the measure <E> which describes the correlations (or fluctuations) of a single particle variable 
x such as the particle energy or transverse momentum. As observed in [22] ], x can also characterize the particle sort. 
Then, x = 1 if the particle is of the sort of interest, say the particle is strange, and x = if the particle is not of this 

dcf _ 

sort, it is a non-strange particle. We define the single-particle variable z — x — x with the overline denoting averaging 
over a single particle inclusive distribution. In the case of the chemical fluctuations, x is the probability (averaged 
over events and particles) that a produced particle is of the sort of interest, say it is strange. One easily observes that 

~z = 0. We now introduce the event variable Z, which is a multiparticle analog of z, defined as Z = ^Zi^LiO^i — ^Oj 
where the summation runs over particles from a given event. By construction (Z) = 0, where (...) represents averaging 
over events. Finally, the <!>— measure is defined in the following way 



We compute $ for the system of particles of two sorts, a and b, e.g. strange and non-strange hadrons. Xi — 1 when 
i— th particle is of the a type and Xi = otherwise. The inclusive average of x and x 2 read 

x = xp x = Pi > a? = x2p * = Pi ' 

x=0,l x=0,l 

where Pi is the probability (averaged over particles and events) that a produced particle is of the a sort. Thus, 

Pl = w 

(N a ) + (N b ) ' 

with N a and Nf, being the numbers of particles a and 6, respectively, in a single event. One immediately finds that 
z = while 

-2_ p p2_ WW (0) 
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where N — N a + N b is the multiplicity of all particles a and & in a single event. 
Since the event variable Z equals N a — xN, one gets 

(Z) = (N a ) - x(N) = , 

{Z 2 ) = {N 2 a ) - 2x{N a N) + x 2 {N 2 ) . 

The latter equation gives 

(Z 2 )(N) 2 = (N b ) 2 (N 2 ) + (N a ) 2 (N 2 ) - 2 (N a )(N b )(N a N b ) , 

which can be rewritten as 

^ = ~ W 2 ) + ^«*?> - (3) 

-2™((«>-WW)- 

The fluctuation measure $ is completely determined by eqs. (|^, |^). So, let us consider its properties within three 
simple models of the multiplicity distribution. 

1) The distributions of particles a and b are poissonian and independent from each other i.e. 

(N 2 ) - (N t ) 2 = {NJ , (4) 
(N a N b ) - (N a )(N b ) , 

where i = a,b. One easily notices that $ — in this case. 

2) The particles a and b are assumed to be correlated in such a way that there are no chemical fluctuations in 
the system. The event chemical composition, which is fully characterized (for a two component system) by the ratio 
N a /N, is assumed to be strictly independent of the event multiplicity. Then, N a = aN and N b = (1 — a) N with a 
being a constant smaller than unity. Since -/V a , N b and N are the integer numbers, a has to be a rational fraction. 
Then, we have 

(N a ) =a(N), (N a ) = (1 - a) (N) , 

<^> - W = a 2 ((N 2 ) - (iV) 2 ) , (N 2 ) - (N b ) 2 = (1 - a) 2 ((N 2 ) - (N) 2 ) , 

(N a N b ) - (N a )(N b ) = o(l - a) ((TV 2 ) - (iV) 2 ) . 

Consequently, (Z 2 ) = and 



(N a )(N b ) r— 

Y (iv) 2 = ~^ ( ~ a) • ( ) 

One sees that the measure is negative (but larger than —1/2) when the chemical fluctuations vanish in the system. 

3) The particles a and b are identified with the positively and, respectively, negatively charged hadrons. Then, the 
charge conservation leads to the strict correlation of the particle numbers: 

N+ - iV_ = Q , 

where Q is the electric charge of the system. In this case we have 

(N+) = (N-) + Q , 

(iV 2 ) - (iV+) 2 - (N 2 ) - (7V_) 2 , 
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(N + N.) - (N+)(N-) = (N 2 _) - (iV_) 2 

Therefore, 

-_ ((M)+Q)(iV_) 
/72\ n2 

«J\£) - (TV-) 2 ) . 



(TV) <iV)3 

When Q = we reproduce the result (^) corresponding to a = 1/2. 

After the illustrative examples let us compute $ for the equilibrium gas which is a mixture of the particles a and 
b. Then, 

(N i ) = X i -^-\nE(V,T,X a ,X b ), (6) 

(AUV b ) - (N a )(N b ) = A a A b in5(V,T, A a ,A b ) , 

OA b OA a 

a ^ 2 



(AT/) - (Ni) = I Ai— 1 lnS(t/,T,A Q ,A b ) , 

where S(V,T, A a ,A b ) is the grand canonical partition function with V, T and Ai denoting, respectively, the system 
volume, temperature and the fugacity which is related to the chemical potential fa as Ai = e"^' with (3 = T _1 . 
When the gas of interest is the mixture of the two ideal quantum gases, the partition function is ]24[] 

lnZ(V,T,\a,\ b )=±9aV J ^\n[l±X a e-^]±g b V J ^ ln[l ± A b e"^] , (7) 

where gi denotes the number of the particle internal degrees of freedom; Ei = \J mf + p 2 is the particle energy with 
nii and p being its mass and momentum; the upper sign is for fermions while the lower one for bosons. 
Substituting the ideal gas partition function (^) into eqs. ([}]), one easilly finds 

< N < } = °> V lw?K^l' (8) 

(N a N b ) = (N a )(N b ) , 

d 3 p ATV* 



(N 2 )-(N i ) 2 = g i V 



(2tt) 3 (Xr 1 e fi E *±l) a 



where, as presviously, the index i labels the particles of the type a or b. It is worth noting that the system volume V 
which enters eqs. (||) cancels out in the final expression of $. Therefore, the measure $ is, as expected, an intensive 
quantity. One observes in eqs. (|8|) that 

(Nf) - (n,) 2 < m 

for fermions, 

(N 2 ) - (N z ) 2 > (Ni) 

for bosons, and 

(N 2 ) - (^> 2 = {N t ) 
in the classical limit where A^ 1 ^> 1. Therefore, one finds from eqs. (B ||) that 

< z 2 and $ < 



(N) 
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when the particles a and b are fermions, 

(Z 2 ) ^ 



(N) 

when the particles a and b are bosons, and 

(Z 2 ) 



(N) 



> z 2 and $ > 



and $ = 



when the particles of both types can be treated as classical. If the particles a and b are of different statistics, the sign 
of <E> is determined by the sign of the expression 

(N a ) 2 ((N 2 ) - (N b ) 2 - (N b )) + (N b ) 2 ((N 2 ) - (N a ) 2 - (N a )) , 

which is either positive or negative depending of the particle masses, their chemical potentials, and the numbers of 
the internal degrees of freedom. 

If all particles are massless and their chemical potentials vanish, the calculations can be performed analytically. In 
this case eqs. (||) give 

5«C(3) f3/4\ T _ 3 „ /0.09\ _ 



where is the Riemann zeta function (C(3) = 1.202); as previously the upper case is for fermions and the lower 
one for bosons. If the particles a and b are both fermions or both bosons, eqs. Q ||) get the form 

— 7 9a9b 



{9a + 9b) 2 ' 

71-2 9a9b /2/3 



(TV) 6C(3) (ga+9b) 2 V 1 

and consequently 



-0.045\ 



0.170y.g Q +. 9b 



Let us now consider the fluctuations in the system of pions and kaons. To be specific, the particles a are identified 
with 7r~ while the particles b with K + or K~ . As we will see, the fluctuations in the n~K + system can be very 
different from those in tt~K~ one. The systems ir + K + and ir + K~, which are not discussed here, are analogous to, 
respectively, tt~ K~ and n~K + . At first we treat the pions and kaons as a mixture of the ideal gases of tt and K. 
Since the pions and kaons are of a given charge (plus or minus) g^ = gx = 1. The masses are taken, respectively, 140 
and 494 MeV. $ as a function of temperature has been computed numerically from the formulas (0J|,^) combined 
with (^). The results, which are obviously the same for the ir~ K + and ir~K~ systems, are shown with the dashed 
lines in Figs. 1-4. The calculations have been performed for several values of the chemical potentials of pions and 
kaons. In the case of pions, fi v = when the system is in the chemical equilibrium. (We obviously neglect here a tiny 
effect of the electric charge conservation.) The finite value of \xk appears even in the equilibrium system of zero net 
strangeness due to the simultaneous baryon and strangeness conservation. For example, the estimated value of fi K 
for the strange (not antistrange) mesons is 38 MeV at T = 160 MeV for the equilibrium hadronic system produced 
in heavy-ion collisions at CERN SPS @. 

In Figs. 1 and 2 we observe a dramatic increase of <& with the temperature. $ also grows with \ik (at fi„ = 0) 
while the dependence on /i^ changes with the temperature. Below T = 120 MeV $ is a decreasing function of but 
above this temperature $ grows with (/i^- is fixed and equals zero). Such a behaviour can be easily understood. 
$ can be approximated as 
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for (N.„) 3> (Nk) which holds for sufficiently low temperatures. (We also assume here that (N^) — (Nk) 2 is not 
larger than (N%) — (N^) 2 .) The growth of \ik leads to the increase of the numerator of the expression while the 
growth of fin enlarges the denominator. At higher temperatures the numbers of pions and kaons are comparable to 
each other and the pion dispersion (N 2 ) — (N^) 2 provides a significant contribution to $. Then, $ grows with /j, n . 

It is a far going idealization to model a fireball at freeze-out as an ideal gas of pions and kaons. A substantial 
fraction of the final state particles come from the hadron resonances. We take them into account in the following way. 
Since the resonances are relatively heavy, their phase-space density is rather low. Consequently, the resonances can 
be treated as classical particles with the poissonian multiplicity distribution Then, one easily shows that 

(Ni) = (Nl) +J2br(N r ) , l = 7T,K (11) 

r 

(N 2 ) - (N t ) 2 = (N' 2 ) - (AO 2 + $>(^> , 

r 

(N„N K ) - (Nk) (Nk) = (KN' K ) - (K){N' K ) + ]T b*(K) , 

r 

where N![ is the number of 'direct' pions or kaons which are described by the formulas (^J); the summation runs over 
the resonances which decay into pions or kaons; N r is the number of resonances of type r and b r is the branching 
ratio of the resonance decay into the pion or kaon channel. The resonances are assumed to decay into no more than 
one pion or one kaon of interest. We denote with * the resonances, such as A"*, which decay into the pion-kaon pair 
under study. 

In the actual calculations we have taken into account the lightest resonances: p(770), lj(782) and AT* (892). Now, 
an important difference between the correlations in the ir~ K~ and tt~ K + system appears. The decays of K*° into 
K + n~ produce the correlation in the tt~ K + system. Analogous correlation in the 7r~ K~ system is absent. $ as a 
function of temperature has again been computed numerically from the formulas (|]j^||) combined with eqs. (||) which 
are now supplemented with eqs. (|TT|). The results are shown with the solid lines in Figs. 1 and 2 for the tt~ K~ system 
and in Fig. 3 and 4 for the tt~ K + one. The calculations have been again performed for several values of the chemical 
potentials of pions and kaons. The chemical potential of p and ui has been taken to be equal to while that of K* 
equals fix- One sees that in the case of ir~K~ correlations the presence of resonances does not change the results 
qualitatively although the value of <& is significantly reduced. The case of ir~ K + is changed dramatically due to the 
resonances. The role of the term corresponding to (N^Nk) — (N- k )(Nk) appears to be so important that $ becomes 
negative for sufficiently large temperatures. It is somewhat surprising that $ from Fig. 4 changes its sign at the 
temperature of about T = 110 MeV which is approximately independent of fix- Such a behaviour can be understood 
as follows. Since (N„) > (Nk) in the domain of the parameter values of interest, we expand the expressions (|J) and 
(||) in powers of (Nk) /(N t ). One observes that the first power terms of (||) and (|^) cancel out each other. Therefore, 
<I> = when the second power terms of (||) and (||) are equal to each other. Taking into account that the kaons are 
approximately classical and consequently $ depends on {Ik roughly as e^ K , one indeed finds that the position of 
$ = is approximately independent of [Ik ■ 

At the end we take an effort to estimate $ from the existing experimental data |25|] which appear to be rather 
scarce. Specifically, we consider the system of K® and negative hadrons produced in pp interactions at the energy 205 
GeV, which is close to the currently available energies of heavy-ion collisions at CERN SPS. This case is expected to 
be similar to the tt~ K + system discussed above. One finds in |2[| that: 

(A_)=2.84, (A 2 .) — (N-) 2 = 3.63 



(N K ) = 0.18 , (N-Nk) - (N_)(N K ) = 0.078 . 

Unfortunately, there are no data on (N K ) — (Nk) 2 which gives a dominant contribution to $. If the multiplicity 
distribution of kaons is poissonian i.e. (A|-) — (Nk) 2 — (Nk), we get $ = —0.004. The dispersion of the negative 
hadron multiplicity distribution is known to follow the so-called Wroblewski formula [^6). Applying the formula to 

the kaons we have (N%) - (N K ) 2 = (0.58(N K ) + 0.29) = 0.16. In this case the kaon multiplicity appears to be 
even narrower than the poissonian one and <& = —0.02. The estimate of $ = 0.007 given in |2^] exceeds the two our 
numbers because the kaon mutiplicity distribution, which is used in |^2| , is (after averaging over the negative hadron 
multiplicity) broader than the poissonian one. We conclude that the existing data give a rather poor information on 
$ in pp collisions. 
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We summarize our study as follows. The measure seems to be a useful tool to study the chemical fluctuations 
in heavy-ion collisions. If the particles of different species are produced independently from each other and the 
multiplicity distributions are poissonian, $ is exactly zero. When the particles are produced in such a way that there 
are no chemical fluctuations (particle ratios are fixed), $ is negative but larger than —1/2. If the nucleus- nucleus 
collision is a simple superposition of N-N interactions the value <& is strictly independent of the collision centrality. 
The same happens when the hadronic system produced in the nucleus-nucleus collisions achieves the equilibrium with 
the temperature and chemical potentials being independent of the impact parameter. The thermal model, which 
seems to be successful in describing the average multiplicities of different particle species, gives a definite prediction 
of <£>, which is positive for bosons and negative for fermions. The correlations in the system of pions and kaons have 
been considered in detail. The estimate of <£> for the it~ K~ system is rather reliable while the prediction concerning 
the ir~K + correlations is sensitive to the details of the model. Since the experimental value of <f> in pp interactions 
can be hardly extracted from the existing data, the fluctuation measurements of nuclear collisions should start with 
the nucleon-nucleon case. 



I am very indebted to ECT* in Trento where the idea to analyze the equilibrium chemical fluctuations was born. 
Numerous fruitful discussions with Marek Gazdzicki, who initiated this study, are also gratefully acknowledged. 
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Figure Captions 



Fig. 1. $— measure of the ir~ K~ correlations as a function of temperature for three values of the pion chemical 
potential. The kaon chemical potential vanishes. The resonances are either neglected (dashed lines) or taken into 
account (solid lines). The most upper dashed line on the right hand side of the figure corresponds to fi v = 100 MeV, 
the central one to fi v = 0, and the lowest line to fj, v = —100 MeV. At sufficiently small temperatures the respective 
dashed and solid lines coincide. 

Fig. 2. measure of the -k~K~ correlations as a function of temperature for three values of the kaon chemical 
potential. The pion chemical potential vanishes. The resonances are either neglected (dashed lines) or taken into 
account (solid lines). The most upper dashed line corresponds to fiK = 100 MeV, the central one to = 0, and the 
lowest line to = —100 MeV. At sufficiently small temperatures the respective dashed and solid lines coincide. 

Fig. 3. The absolute value of $— measure of the n~K + correlations as a function of temperature for three values of 
the pion chemical potential. The kaon chemical potential vanishes. The resonances are either neglected (dashed lines) 
or taken into account (solid lines). The most upper dashed line on the right hand side of the figure corresponds to 
yUjr = 100 MeV, the central one to = 0, and the lowest line to ^ = —100 MeV. At sufficiently small temperatures 
the respective dashed and solid lines coincide. 

Fig. 4. The absolute value of <&— measure of the n~ K + correlations as a function of temperature for three values 
of the kaon chemical potential. The pion chemical potential vanishes. The resonances are either neglected (dashed 
lines) or taken into account (solid lines). The most upper dashed and solid lines correspond to [Ik = 100 MeV, the 
central ones to [ik = 0, and the lowest lines to /j,k = —100 MeV. At sufficiently small temperatures the respective 
dashed and solid lines coincide. 
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